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We show that QCD Green functions in Landau Gauge exhibit sizable 0(l//i^) corrections to the expected 
perturbative behavior at energies as high as 10 GeV. We argue that these are due to a (A^)-condensate which 
does not vanish in Landau gauge. 



1. Computing Green functions 

1.1. The lattice settings 

We have performed a series of run in pure Yang- 
Mills SU3 with [3 = 6.0, 6.2, 6.4, 6.8 i.e. a'^ rang- 
ing from 0.03 to 0.1 fm. The lattices were hyper- 
cubic: 16^, 24"^ and once 32^. The statistics was 
of 1000 configurations per run. More details are 
available in H,^. 

1.2. The method 

The method used to compute the strong cou- 
pling constant is described in It consists 
in computing non-perturbatively the renormal- 
ized Green functions in MOM schemes. For a 
given scale /z the gluon fields are renormalized so 
that the gluon propagator at that scale is equal 
to l//i^ (times the standard tensor in Landau 
gauge). The renormalized coupling constant is 
defined so that the three-gluon vertex function. 



projected on the tree-level tensor, at that scale, 
is equal to the renormalized coupling constant. 

In fact two different MOM schemes have been 
used. In the "symmetric" one, the three-gluon 
vertex function is taken with gluon momenta = 
P2 = pi = A*^: while the "asymmetric" one, called 
MOM corresponds to Pi = p| = H^,P3 = 

1.3. Results without power corrections 

The three-gluon vertex function as defined 
above gives directly a nonperturbative value for 
as in) for all values of /i such that /i^ — J2^Pt- 
The gluon propagator evolves at large momenta 
according to the perturbative formula. This evo- 
lution depends of course on the value of as (n) at 
some properly chosen scale fj, (large enough to be 
in the perturbative regime) in the range covered 
by the lattice calculation 0. In this subsection 
we neglect all power corrections. Thus we as- 
sume that, for large enough /i the evolution of 
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the propagator is fully described by the pertur- 
bative formula taken to three loops. Once as is 
known we compute Aqcd to three loops in the 
considered scheme, and we eventually translate it 
to the standard Aj^. Of course, if we really are 
in the perturbative regime A^|g must be constant 
as the energy varies. 

We obtain from the three-point Green func- 
tions 



a^°^(4.3 GeV) 



0.269(3); Ams(4.3 GeV) = 299(7) 



^MOM^g g ggy^ ^ 0.176(2); Aj^(9.6 GeV) = 266(7) 

and from the two point-Green function 

^MOM^g g ggy^ ^ 0.193(3); A5^(9.6 GeV) = 319(14) 

We see that there is a problem in that Aj^ 
varies much more than allowed by the statistical 
errors. 

2. Power corrections. 

We now add in our fits power corrections 
to the three loop formula. In Q the strong cou- 
pling constant is fitted according to 



asifi^) = Q;s(p^) Sloops + 



(1) 



in a large energy window: 2 to 10 GeV. The result 
of the fit is 



Aj^ = 237 ± 4 ± lOMeV; c : 



0.63 T O.OSGeV^ (2) 



in strikingly good agreement with the result from 
the ALPHA collaboration 1^ who has used a to- 
tally different method. The effect of power cor- 
rections on the expected constancy of Aj^ is il- 
lustrated in fig |l| 

2.1. Povirer corrections and OPE 

What can be the theoretical reason for these 
large power corrections ? In the theoretical 
framework of OPE a quadratic power requires a 
condensate of a dimension- two operator. Since 
we are working in the Landau gauge, such an op- 
erator exists: < A^A^ >. Interestingly enough it 
is the only dimension two candidate for a vacuum 
expectation value in the Landau gauge. 




• A(a), p=6.0 (16*) 
A(a), p=6.2 (16*) 

• A(a), p=6.4 (16*) 
A(a), |i=6.8 (16*) 
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Figure 1. The upper plot shows Aj^ computed 
without power corrections while the lower plot 
shows Aj^ computed form the perturbative con- 
tribution to as once the power correction has 
been substracted, eq. (|]) . 



We start from the OPE expansion of both the 
gluon propagator |6j and the symmetric three- 
gluon Green function, keeping only the relevant 
terms: 



T 



(3) 



(c2)^ta'&'(P) 



Al,{0) A^(0) 



T(Alip,)AiiP2)A'i{ps) 



(rfo)X(Pl'P2,P3)' 



+ id2C%rUpi,P2,P3) : <.(0) A'UO) :+■■■ 

The Wilson coefficients depend on the renor- 
malization scale. The leading coefficients, co{n) 
and do(/^) are exactly the perturbative estimate 
of the considered Green functions which is known 
to three loops. The first subleading coefficients 
C2(/^) and d2{p) are proportional to 1//^^ for di- 
mensional reasons and they also depend logarith- 
mically on the renormalization scale. We have 
computed to the leading logarithms approxima- 
tion the anomalous dimensions of these coeffi- 
cients and that of the A^ operator. 

We have then performed the following test of 
OPE: A simultaneous fit was first done of the two- 
point Green function and of Ug deduced from the 
symmetric three-point Green functions computed 



(4) 
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on the lattice. We have used the known depen- 
dence of the Wilson coefficient on the scale but, at 
this stage, fitted two different values for < >. 
We then check the consistency of the resulting 
values of < > . 



the asymmetric as deduced from the asymmetric 
three-gluon Green function leads to: 

{A^)prop - (1.39(14)GeV)' ; {A^)aipha - (2.3(6)GeV)' 



Ipha 



{A')prop - (1.55(17)GeV)^; M^)a/p/.a = (1.9(3)GeV)^ ^(A^ 



= 1.7(3); Aj^= 260{18)MeV 



Ipha 



= 1.21(18); Aj^ = 233(28)MeV 



The two fitted values of < > agree within 
one standard deviation. This is a quite interesting 
result which tends to confirm the hypothesis that 
the A'^ condensate is at the origin of the power 
corrections. It is interesting to stress that if we 
use only the two loops formulae for the leading 
coefficients the fitted condensates differ by more 
than twenty standard deviations. This is a les- 
son for the use of OPE: if the leading coefficients 
are not expanded far enough in perturbation|^ the 
OPE cannot be applied. 

2.2. OPE for the asymmetric three-gluon 
Green function 

OPE is more tricky to apply to the asymmetric 
Green function because of the presence of a zero 
momentum gluon. To be short, this mixes up the 
distinction between hard and soft gluons which is 
at the basis of the OPE approach. The issue has 
been discussed in Strictly speaking OPE ap- 
plied to this operator introduces other operators 
than the identity and the considered above. 
One needs to consider gluon to vacuum matrix el- 
ements of operators with an odd number of gluon 
fields. 

In order to be able to relate nevertheless the 
asymmetric three-gluon Green function measured 
on the lattice to the condensate the authors 
of have introduced a "vacuum insertion hy- 
pothesis" (or factorization hypothesis) consisting 
in taking the matrix elements of three-gluon op- 
erators as the products of matrix elements of 
one gluon operators and of two-gluon operators. 
Then, the only non trivial matrix element is again 
< > and one can perform the same exercise as 
above. The combined fit of the propagator and 



The agreement between both estimates of 
< > is less satisfactory than in the previous 
case, maybe due to the added assumption of fac- 
torization or to the necessity of going to four loops 
for the dominant coefficients. 

As a conclusion, we would like to emphasize 
the role of lattices as a beautiful tool to study 
the QCD vacuum and in particular the vacuum 
condensates and also to test the practical applica- 
bility of OPE. Notice also [|j that there is a good 
evidence for power corrections in the gluon Green 
functions computed with dynamical quarks. 
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